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derived,  followed  by  the  infinite  cylinder  confutation  which  should  be 
accurate  for  wavelengths  short  compared  with  cylinder  length.  In  order 
that  the  computations  may  be  extended  into  the  infrared  and  visible 
regimes,  it  is  necessary  to  incorporate  the  optical  properties  of  the 
fibers.  ,  We  do  this  by  employing  the  Drude  Model  for  conductivity  (or 
complex  /dielectric  constant) ,  and  also  introducing  the  dependence  of 
conductivity  on  both  fiber  diameter  and  electron  mean  free  path. 
Numerical  results  are  presented  for  copper  and  lead  fibers ,  and  the 
approximate  methods  are  seen  to  agree  very  well  with  the  variational 
computation  in  those  limiting  regions  where  they  should  apply. 

^ Curved  fibers  are  considered  in  this  report.  An  exact  integral 
equation  is  derived  for  the  general  case,  and  sane  approximate  results 
are  then  given  for  special  fibers.  Also  discussed  are  the  optimum 
conditions  for  target  obscuration  by  a  cloud  of  particles  from  a  mass 
efficiency  standpoint.  Particle  optimization  studies  are  performed 
under  varying  requirements  on  absorption,  reflection  and  transparency 
of  the  particle  cloud. 
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1.0  INTRODUCTION 


During  the  course  of  the  past  several  years,  a  theory  based 
on  the  variational  method,  along  with  associated  computer  codes, 
has  been  developed  at  Panametrics  for  analyzing  the  electromagnetic 
scattering  and  absorption  from  thin  conductive  fibers  of  arbitrary 
size,  conductivity  and  orientation.^  Extensions  and  refinements 
of  this  theory  have  now  been  completed  and  programmed,  as  described 
herein. 

We  begin  by  summarizing  the  basic  equations  used  in  the 
variational  computation  for  arbitrary  fibers.  The  quasistatic 
model  appropriate  at  long  wavelengths  is  then  derived,  followed 
by  the  infinite  cylinder  computation  which  should  be  accurate  for 
wavelengths  short  compared  with  cylinder  length. 

In  order  that  the  computations  may  be  extended  into  the 
infrared  and  visible  regimes,  it  is  necessary  to  incorporate  the 
optical  properties  of  the  fibers.  We  do  this  by  employing  the 
Drude  model  for  conductivity  (or  complex  dielectric  constant),  and 
also  introducing  the  dependence  of  conductivity  on  both  fiber 
diameter  and  electron  mean  free  path. 

Numerical  results  are  then  presented  for  copper  and  lead 
fibers,  and  the  approximate  methods  are  seen  to  agree  very  well 
with  the  variational  computation  in  those  limiting  regions  where 
they  should  apply. 


Curved  fibers  are  considered  in  Section  4.  An  exact  integral 
equation  is  derived  for  the  general  case,  and  some  approximate 
results  are  then  given  for  special  fibers.  In  Section  5,  we 
consider  the  optimum  conditions  for  target  obscuration  by  a  cloud 
of  particles  from  a  mass  efficiency  standpoint.  Finally,  Section 
6  deals  with  particle  optimization  studies  under  varying  require¬ 
ments  on  absorption,  reflection  and  transparency  of  the  particle 
cloud. 

2.0  THEORY 

2.1  The  Variational  Equations 

The  variational  method  has  been  discussed  often*l-4  we  will 

» 

simply  summarize  the  important  equations  here.  The  basic  integral 
equation  is 


with  surface  impedance  Z  given  by 

g  I  (  g  a  ) 

7  =  -  °  —  (  ?  ) 

2nkca  (e”  +  ie')  Ii(ga) 

2  2 

=  k0  [(cos  0O  -  e’)  +  Ie"! 
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For  the  constant  A  appearing  in  Eq.  (3),  the  variational  method 
then  gives 
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all  of  which  can  be  evaluated  explicitly. 


The  extinction,  absorption  and  scattering  cr o s s- s e c t i ons  are 


then  given  explicitly  by 
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2.2  The  Quasistatic  Model 


In  the  following,  we  re-derive  the  equations  of  Ref.  5  and  6 
for  absorption,  scattering,  and  extinction  worked  out  in  1965. 
Although  the  previous  equations  are  satisfactory  at  radar  freque¬ 
ncies,  one  must  consider  the  relaxation  effects  of  Sections 
2.4  and  2.5  when  calculating  optical  and  infrared  electromagnetic 


interactions. 


M 


First  consider  the  electric  field  E^  inside  a  long,  thin 
spheroid,  given  by 


Ei  =  E0  -  L  4n  P. 


where  L  =  depolarizing  factor,  P  =  polarization,  and  E0  =  homo¬ 
geneous  applied  field.  In  the  case  of  a  long,  thin  spheroid, 
which  is  a  good  approximation  to  our  conductive  filament,  the 
depolarizing  factor  is 


L  =  4(f) 


where  a  =  particle  radius  and  i -  particle  length. 


The  definition  of  the  polarization  is 


( e -1 )  E. 

1 


where  e  =  normalized  dielectric  constant.  Combining  (9)  and  (10) 
we  see  that  the  internal  field  is  given  by 


Ei  =  — 


L  (e-i: 


We  next  take  e  to  be  complex,  i.e.,  e  =  e'  +  ie  .  The  internal 
field  then  becomes 


Ei  =  E0 


I  1  +  L  (  e  '  - 1 )  Le  -  i 
<  [L(e'-l)  +  l]2  +  [Le"] 


-  iLs  -  iLs 
2  .  ,,,2 


and  we  see  that,  in  the  quasistatic  approximation,  we  have  both 
an  in-phase  and  a  quadrature  component  of  the  internal  field. 


TOW 


5 


The  electric  dipole  moment  of  our  Rayleigh  particle  is 


defined  in  terms  of  the  polarization  as 


(1)  =  py, 


where  Vp  =  particle  volume.  The  polarizability  a  of  the  particle 
is  related  (by  definition)  to  the  dipole  moment  as  follows: 


(1)  =  a  E0  . 


Combining  (12),  (13),  and  (14),  we  obtain  the  needed  expression 

for  the  complex  electric  polarizability  of  the  particle: 
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The  well  known  extinction  theorem**  relates  the  extinction 
cross  section  of  any  scatterer  to  its  normalized  forward  scattering 
amplitude,  S(0)  by  the  following  relation: 


®ext  =  — o  Re  [S ( 0) }  , 


where  oext  =  extinction  cross  section,  k  =  2 n/k0,  X0  =  free  space 
wavelength  of  the  incident  wave.  Correct  to  order  k*>  ,  S(0)  is 
related  to  the  polarizability  by** 


S( 0)  =  ik3  a  +  3  k6  a2 , 


where  the  complex  polarizability  is  given  by  Eq.  (15). 
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The  absorption  and  scattering  cross  sections  are,  respec 
tively,  given  by 


°abs  =  4  it  k  Re  (-ia) 


and  csca  =  j  nk4  I  a.  I  2  , 


where  a  is  given  by  Eq.  (IS).  After  considerable  manipulation 
the  following  expressions  for  the  cross  sections  emerge: 


ext  =  kV  p  [  B  +  (B2  +  A2)] 


:abs  =  kVpB 


-  <a2  +  B2> 
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and,  in  summary  from  the  preceding  sections  plus  Sec.  2.4  and  2.5 


e'  =  1  -  21 
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o  =  aG  11  -  |  (1-6)  (£)]  (28) 

In  an  attempt  to  analyze  the  effects  of  Eq s .  (25)  and  (26) 

on  the  absorption  cross  section,  we  have  substituted  these 

expressions  into  Eq.  (24),  and  then  inserted  Eq.  (24)  into 
Eq .  (21).  The  resulting  expression  for  absorption  cross  section 

is  too  complicated  for  easy  analysis.  However,  there  is  no  doubt 
that  in  the  regime  of  (  w  x  )  2  <  <  1  (  the  absorption  cross  section 

is,  in  general,  diminished  from  its  value  without  the  inclusion 
of  relaxation  effects.  It  turns  out  that,  for  copper,  wt  =  1 

corresponds  to  X0  ±  30  microns. 

We  have  recently  shown  that  (at  least,  for  the  filament 
types  of  interest  here),  we  can  greatly  extend  the  region  of 
applicability  of  the  q ua s i s t a t i c  theory  by  letting  £=  n/k0  in  the 
k  0  f  >  1  regime.  We  do  this  only  in  the  depolarizing  factor,  L', 

given  by 


L'  =  4 


k  a 

o 


[  £  n  (r2— )  -1] 
k0a 


(29) 


Note  that  this  substitution  cannot  be  used  to  predict  the  scat¬ 
tering  cross  section  when  k0£ >  n ,  and  is  good  only  for  absorption. 
We  will  comment  shortly  on  the  use  of  this  formula. 


2.3  The  Infinite  Cylinder  Appr o xima t i on 

Using  the  results  of  Wait7,8  for  scattering  by  infinitely 
long  cylinders  of  arbitrary  dielectric  constant,  we  have  derived 
the  expression  for  absorption  cross  section  per  unit  length  of 
the  infinite  cylinder.  The  general  expression,  obtained  by 
considering  the  inward  radial  component  of  the  Poynting  vector  on 
the  lateral  surface  of  the  cylinder,  is 


VJio  ,  k 

eQ  sin2  Re  [(-^)  J  (k*a)J  ’(ka)] 


.  y  k  a  s  i  n  © 

J0(ka)  _  (~)  sin  ©D  [  £n  ( - j - -)  +  y  +  j— ]  ( ka  )  J0  '  (  ka  ) 

k2  =  ne  kc2 


Y  =  .5772 


We  have  compared  the  results  of  our  high  frequency  asymptotic 
quasistatic  expression  for  the  absorption  cross  section  with  this 
expression.  The  results  in  all  (six  or  seven)  cases  were  in 
close  ( ~  20%)  agreement,  as  will  be  shown  below.  This  close 
agreement  was  both  surprising  and  gratifying  to  us.  It  is  now 
evident  that  we  have  a  quasistatic  theory  which  can  predict  the 
absorption  cross  section  of  thin  conductive  filaments  over  a  very 
wide  range  of  k£,  and  which  predicts  the  scattering  cross  section 
in  the  region  k0£  <  i. 


2.4  The  Drude  Model 


Our  early  workS,  6  the  area  of  target  obscuration  was 

restricted  to  the  microwave  range  of  incident  electromagnetic 


energy.  We  did  not  at  that  time  anticipate  that  the  infrared  and 
visible  regions  of  the  spectrum  would  be  applicable  to  the 
absorbing  particle  cloud  techniques  which  we  developed.  Although, 
during  our  theoretical  programs  with  the  U.S.  Army  Chemical 
Systems  Lab,  such  applications  were  being  explored,  we  were 
so  involved  in  developing  a  comprehensive  theory  that  we  did  not 
include  the  opt i c a 1  properties  of  the  (metallic)  filaments  whose 
cross  sections  we  were  calculating.  Indeed,  it  has  not  been 
until  the  present  AFOSR  program  that  we  have  included  these 
properties.  In  this  and  the  following  sub-section  we  deal  with 

(1)  the  inclusion  of  the  optical  properties  in  both  theories,  and 

(2)  the  inclusion  of  the  dependence  of  electrical  conductivity 
upon  particle  radius  and  the  electron  me a n- f r e e-pa t h. 

In  a  recent  paper  by  Ordal,  et  a 1^,  the  application  of  the 
Drude  model  for  the  prediction  of  complex  optical  dielectric 
constant  was  compared  with  measured  values  of  the  real  and 
imaginary  parts  of  the  optical  dielectric  constant  for  a  number 
of  metals  (Al,  Cu,  Au,  Fb,  Ag,  and  W)  .  Tabular  experimental 
results  are  also  given  for  Fe,  Pt,  Co,  Ni,  Ti  and  Pd.  This 
model,  which  is  based  on  the  free  electron  theory  of  metals,  is 
in  surprisingly  good  agreement  with  the  observed  experimental 
results.  We  realize,  that,  for  certain  transition  elements  such 
as  Fe,  the  model  has  drawbacks.  For  such  cases,  one  must  resort 
to  the  use  of  tabular  experimental  data.  For  our  present  discus¬ 
sion,  we  will  choose  Cu  as  the  substance  comprising  our  fibers, 
and  will  utilize  the  Drude  model  in  the  calculation  of  the 
various  electromagnetic  cross  sections.  A  good  exposition  of 


this  is  given  in  Wooten's  booklO,  in  which  the  normalized  complex 
dielectric  constant  is  derived  (note  that  Eq.  (32)  corrects  a 
mistake  on  p.  53  of  Wooten's  book): 


e  =  e  '  +  i  e 


(31) 


u>  2  x  2 

e'  =  i - 2 - x  (32) 
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_ E _ 
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Dsing  MKS  units,  the  plasma  frequency  Up  is  given  by 


(i) 


P 


2  = 


(34) 


in  which  N  =  electron  density  (m_3),  e  =  electronic  change,  m  = 
effective  mass  of  the  electron,  and  e0  =  permittivity  of  free 
space  =  ( 1  /  3  6  tt  )  x  10~9  farads/m. 

The  quantity  t  is  the  electron  relaxation  time,  which  is  the 
time  required  for  randomization  of  the  momentum  vector  of  an 
electron  in  the  (metallic)  lattice.  For  our  purposes,  it  is 
instructional  to  cast  the  dielectric  constant  in  terms  of  the 


low  frequency  electrical  conductivity,  a,  given  by 


From  (34)  and  (35),  we  obtain 


Using  Eq.  (36)  in  (31)  and  (32),  we  have 


e  * 


1 


+ 


Ux)2 


(36) 


(37) 


to  e , 


[1  +  (  <i)T  )  ] 


(38) 


Typical  values  of  the  relaxation  time  r  are  on  the  order  of 
10-14  sec. 


It  is  easily  shown  that  Eq.  (37)  can  be  written  in  terms  of 


e '  =  1  -  e"(ux)  .  (39) 

Thi s  equation  show  s  us  that,  for  all  frequencies  significantly 
below  the  visible  and  infrared  (i.e.  the  microwave  region),  le'l 
<<  le"l.  And,  from  Eq.  (38),  we  see  that  e"  goes  to  its  low 
frequency  value  e"  =  o/we0  for  (wt)2  <<  1. 

The  reason  for  the  above  analysis  is  to  determine  whether  or 

not  the  Drude  model  can  be  utilized  at  1 ow  frequencies.  Although 

the  low  frequency  asymptotic  value  of  e'  differs  significantly 

from  a  value  of  unity,  which  is  normally  assumed  for  metals  at  1 ow 

frequency,  the  ratio  le"/e'l  will  always  be  very  large  when  («dt)2 
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<<  1.  Therefore,  the  use  of  the  Drude  model  throughout  the 
region  10~6  m  <_  £.  10“^  m  appears  to  be  justified,  and  we  feel 
confident  in  using  Eqs.  (37)  and  (38)  in  the  derivation  of  the 
electromagnetic  cross  sections  throughout  this  entire  wavelength 
range  . 

2.5  The  Reduced  Conductivity 

When  one  or  more  dimensions  of  a  conductive  material  (metal 
or  semiconductor)  are  on  the  order  of  the  mean  free  path  of  the 
conduction  electrons,  electron  collisions  with  the  surface  will 
significantly  reduce  the  mean  free  time,  and  hence  the  mean  free 
path  A  given  by 

A  =  vpT  ( 40 ) 

where  vp  =  Fermi  velocity  and  x  =  relaxation  time  discussed 
previously.  Since  the  electrical  conductivity  is  proportional  to 
x  (Eq.  35),  the  conductivity  will  be  reduced. 

The  classic  work  on  this  subject  was  done  in  1938  by  Fuchs. H 
In  a  more  recent  paper,  Dinglel2  reviews  the  subject  and  provides 
some  useful  numerical  computations.  The  key  equation  in  Dingle's 
paper  is  his  Eq.  (2.3)  which  relates  the  effective  conductivity  a 
to  the  bulk  conductivity  o0,  as  a  function  of  the  mean  free  path 
A,  the  wire  radius,  a,  and  the  quantity  s  that  is  the  probability 
of  an  elastic  collision  at  the  surface: 


A  value  of  e 


1/2  is  frequently  used  as  the  surface  scat¬ 


tering  coefficient.  Taking  this,  we  have 


°  U  ~  16  (i)]  *  (42) 

We  have  had  difficulty  in  obtaining  numerical  values  for  mean 
free  paths  from  the  literature.  However,  Kittell3  provides  a 
good  background  as  well  as  quantitative  data  for  a  number  of 
metals.  Taking  copper  as  the  subject  material  a  value  of  A  *  4.2 
x  10_8m  is  given  in  Table  10.1  of  Rittel’s  book.  Using  Eq.  (40), 
and  taking  the  Fermi  velocity  Vp  =  1.6  x  10^  m/sec  (Kittel, 
p.  240),  we  obtain  a  mean  free  time  of  x  =  2.6  x  10~14  sec.  This 
is  in  fairly  good  agreement  with  the  value  x  =  1.9  x  10~14  Sec, 
which  we  obtain  from  Ref.  9  for  copper. 


Utilizing  the  above  value,  we  obtain  for  copper 


o  =  °0  [1  ~ 


7.9  x  10-9 


(43  ) 


where 

the  radius 

a  has  the 

unit 

of  meters. 

This  equation  shows 

that. 

if  a  =  1.6 

x  1 0  _  8  m 

(160 

An  g s  t  r  om  s  )  , 

the  conductivity  is 

roughly  half  its  bulk  value.  If  the  radius  is  0.1  micron,  the 
conductivity  is  92%  of  the  bulk  value. 


This  exercise  was  done  to  show  that,  indeed,  one  must 


consider  the  particle  size  effect  upon  electrical  conductivity. 


9 
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when  computing  the  absorptive  and  extinction  properties  of  thin 
metallic  fibers  or  films. 

3.0  RESULTS  FOR  COPPER  AND  LEAD  FIBERS 

All  of  the  above  work  has  been  reduced  to  operating  computer 
codes  with  quantitative  graphical  outputs.  It  should  be  noted 
that,  as  a  part  of  our  present  AFOSR  contract,  we  have  received  a 
Hewlett  Packard  HP9000  Model  520  computer.  This  machine  will 
replace  our  present  HP9835  computer,  which  presently  requires  50 
hours  to  complete  a  full  set  of  scattering,  absorption,  and 
extinction  data  on  a  given  particle  over  a  wavelength  range  from 
10  cm  to  1  micron.  The  new  machine  will  reduce  this  time  by  a 
factor  of  between  100X  and  1000X. 


As  a  first  example,  consider  copper.  We  take  fiber  radius 
of  5  x  10“8  m,  unreduced  bulk  conductivity  of  5.8  x  10?  mho/m, 
relaxation  time  1.9  x  10"^^  sec,  mean  free  path  of  3  x  10-^  m. 
The  reduced  conductivity  can  now  be  computed  from  Eq.  (42),  and 
using  that  result  the  complex  optical  dielectric  constant  is 
obtained  from  Eq s .  (37)  and  (38).  Results  of  this  computation 
are  plotted  in  Fig.  1  vs.  wavelength  over  the  range  from  10 
cm  to  1  mi  cron . 

Knowing  e'  and  e  ",  the  surface  impedance  of  the  fiber  may  be 
computed  from  Eq.  (2).  The  results,  which  are  again  frequency- 
dependent,  are  shown  in  Fig.  2. 


1  5 


Choosing  a  fiber  length  of  5  x  10~4  m,  the  extinction, 
absorption  and  scattering  cross-sections  may  now  be  obtained  from 
Eq.  (7)  of  the  variational  method,  and  are  plotted  vs.  wavelength 
in  Fig.  3.  Note  that  both  the  scattering  and  absorption  cross- 
sections  rise  to  peak  values  a  t  A  ~  10-3  then  fall  off.  The 
extinction  cross-section  must  equal  the  sum  of  the  scattering  and 
absorption  cross-sections,  and  this  is  seen  to  be  accurately 
obeyed  except  at  the  shorter  wavelengths  X  10-4  m,  where  we 
believe  that  crs  is  in  error  due  to  employing  too  few  points  in 
the  numerical  integration  (which  affects  only  os). 


In  Fig.  4  we  show  the  cross-sections  obtained  using  the 
quasistatic  theory,  for  the  same  example.  As  expected,  both  ca 
and  os  are  in  excellent  agreement  with  the  previous  values  at  the 
longer  wavelengths.  The  absorption  cr o s s- s e c t i o n  moreover,  is 
seen  to  agree  with  that  of  Fig.  3  for  the  shorter  wavelengths 


also,  providing  good  confirmation  of  the  high 
rizing  factor  proposed  in  Eq.  (29). 


depola- 


In  Fig.  5  oa  is  plotted  using  the  infinite  cylinder  approxi¬ 
mation  Eq.  (30)  appropriate  for  short  wavelengths.  Almost  exact 
agreement  is  found  with  both  preceding  results,  for  X  i.  10~3  m> 
and  hence  for  k££  3,  about  as  one  would  expect. 


Figures  6  and  7  give  respectively  the  variational  and 
quasistatic  results  for  the  same  cylinder,  but  shorter  by  a 
factor  of  ten.  The  peaking  and  crossover  behavior  from  long  to 
short  wavelength  is  now  seen  to  occur  at  X  ~  10"4  m  as  one  would 


anticipate.  Agreement  between  the  two  figures  is  analogous  to 
that  of  the  first  example,  although  the  peak  height  predicted  by 
the  quasistatic  theory  for  oa  is  seen  to  be  somewhat  low  this 
time.  For  the  infinite  cylinder  approximation,  the  curve  of 
Fig.  5  applies  without  change,  and  again  oa  is  in  excellent 
agreement  with  the  variational  values  at  the  shorter  wavelengths 
f  or  which  k H  2  3  . 

Next,  consider  lead  fibers.  The  resulting  dielectric 
constant  and  impedance  are  shown  vs.  wavelength  in  Figs.  8  and  9, 
with  fiber  parameters  as  listed  in  the  figures.  The  variational 
results,  along  with  the  long  and  short  wavelength  approximations, 
are  shown  in  Figs.  10-12,  respectively.  Agreement  is  analogous 
to  that  obtained  with  copper.  Note  that  in  these  cases  oa  shows 
saturation  effects  for  wavelengths  2  10~5  m. 

Figures  13  and  14  show  that  upon  increasing  the  fiber  radius 
by  a  factor  of  10,  the  dielectric  constant  is  substantially 
unchanged,  while  the  surface  impedance  decreases  at  all  wavelengths 
by  about  two  orders  of  magnitude.  The  resulting  cross-sections, 
given  in  Figs.  15-17,  again  show  the  same  pattern  of  agreement 
between  the  variational  (Fig.  15)  and  approximate  results. 

Finally,  Figs.  18  and  19  (which  should  be  compared  with  the 
short  wavelength  approximation  of  Fig.  17)  give  the  corresponding 
values  when  the  fiber  length  is  decreased  by  a  factor  of  ten. 
This  time  the  absorption  cross-section  peaks  at  k  j?  ~  3,  precisely 
as  observed  originally  with  copper. 


The  agreement  seen  consistently  throughout  these  computations 
provides  striking  confirmation  of  the  validity  of  the  variational 
results,  and  also  provides  some  explicit  bounds  on  the  regimes  in 
which  either  of  the  approximate  methods  might  safely  be  used 
independently. 

4.0  SCATTERING  BY  CURVED  CONDUCTIVE  FIBERS 
4.1  Preliminary  Survey 

All  the  numerical  results  obtained  by  Panametrics  to  the 
present  time  for  scattering  and  absorption  by  conductive  fibers 
have  dealt  with  straight  fibers.  The  question  naturally  arises 
then  as  to  the  effects  of  curvature  of  the  fibers  on  the  scattering 
and  absorption  efficiency. 

A  search  of  the  literature  reveals  very  little  work  on 
curved  wires*  all  of  that,  with  two  exceptions,  involving  perfectly 
conducting  wires.  The  curved,  perfectly  conducting  wire  was 
apparently  first  considered  by  Aharoni  in  1946.14  His  equations 
were  applied  to  circular  loop  and  spiral  antennas  by  M  e i . 1 5  jn 
1956  Eouyoumjian  considered  back-scattering  from  perfectly 
conducting  circular  loops. The  two  exceptions  to  the  perfectly 
conducting  case  are  the  work  of  Philipson,  who  considered  lossless 
dielectric  rings,  17  and  Acquista,  who  considered  wavy  cylinders. 

In  both  of  these  latter  cases,  however,  the  scatterer  was  taken 
to  be  only  a  perturbation  on  its  surroundings,  so  that  a  full 
integral  equation  approach  was  not  required. 


We  have  derived  the  integral  equation  for  curved  fibers, 
having  finite  conductivity,  from  first  principles.  The  usual 
thin-iire  analysis  invariably  assumes  that  the  electric  field  can 
be  expressed  in  terms  of  a  current  filament  concentrated  on  the 
axis  of  the  fiber,  and  this  is  physically  somewhat  unsatisfactory 
especially  when,  as  in  the  present  case,  we  must  work  with  both 
electric  and  magnetic  fields  at  the  surface  of  the  fiber.  We  use 
instead  an  approach  based  on  Huygen's  principle, 19  which  states 
that  fields  generated  by  the  tangential  components  of  E  and  H 
(distributed  along  the  surface  of  the  fiber)  must  precisely 
cancel  the  axial  components  of  the  incident  electric  field  along 
the  fiber  axis. 

This  results  in  an  equation  involving  integrals  of  the  two 
unknown  functions  E  and  B  along  the  fiber.  Taking  the  thin-wire 
limit  where  fiber  radius  is  very  small  compared  to  incident 
wavelength,  the  second  of  these  integrals  is  fairly  straight¬ 
forward,  and  for  good  conductors  is  in  t  e  r  pr  e  t  a  bl  e  as  the  field 
due  to  a  distribution  of  surface  currents.  The  first  integral 
behaves  differently,  however.  The  kernel  reduces  to  a  delta- 
function,  resulting  in  a  term  in  E  at  the  field  point  of  evaluation 
of  the  integral  equation.  Surface  values  of  E  and  H  can  then  be 
related  by  a  surface  impedance  concept  to  finally  give  a  pure 
integral  equation  for  the  current. 


4.2  The  Integral  Equation 

When  an  incident  electromagnetic  field  Eiuc  illuminates  a  body 
in  free  space,  the  resulting  fields  are  related  by  Huygens’ 
principle,  which  states  rigorously  thatl9 

Einc(r_)  -  (1/4*)  V  x  jdo'kn'  x  E+ (  r  '  )  g(kR) 

/E(r)  outside  (44) 

-(1/4*)  V*  Vi  J  do '  in '  x  H+ (  r_' )  g(kR)  =  \ 

/  0  inside  . 

The  left-hand  side  (LHS)  of  this  equation  consists  of  the  sum  of 
the  electric  fields  due  to  the  incident  wave  and  surface  distribu¬ 
tions  over  the  body  of  magnetic  and  electric  dipoles,  respectively. 
Here 


g(kR)  =  ( 1/ kR ) e i RR ,  R  =  j  r  -  r '  |  , 


(45  ) 


where  r_*  and  r  are  the  source  point  and  field  point,  respectively. 

Equation  (44)  states  that  the  E  field  is  given  by  the  LHS 
for  all  field  points  outside  the  surface.  On  the  other  hand,  for 
all  field  points  within  the  surface  the  LHS  vanishes  identically, 
i.e.  the  surface  field  distributions  must  precisely  cancel  the 
incident  wave.  This  latter  statement  is  sometimes  known  as  the 
extinction  theorem,  or  the  extended  boundary  condition. 

The  extinction  theorem  is  applied  to  the  curved  fiber,  shown 


in  Fig.  20,  as  follows. 


Let  a 


fiber  radius,  and  p(s) 


radius 


of  curvature  as  a  function  of  position  s  along  the  axis  of  the 


fiber.  We  assume  that 


a  / Pnin  (  ^  1 


(46  a) 


a/b  <  <  1 


(46b) 


ka  <  <  1  .  (46c) 

i.e.,  the  fiber  r v.  j i  u  s  is  much  less  than  the  minimum  radius  of 
curvature,  the  fiber  half-length  b,  and  free  space  wavelength  k  = 
2 n /X,  respectively.  We  also  assume  the  fiber  to  have  moderate 
to  large  conductivity,  so  that  axial  currents  will  be  induced  and 
guided  along  the  fiber. 

Now  requiring  that  the  aiia)  com  pone  nt  of  Eq .  (44)  vanish 

along  the  fiber  axis  gives 

( 1/  4n  )  s  •  v  x  v  x  J"  do '  in'  x  H+  (  r_’  )  g(kR) 

+  <l/4n)  s  .  vx  Jdo'  kn'  x  E+ ( )  g(kR)  =  s  •  E*nc(5)  .  (47) 

Note  that  this  equation  is  still  exact,  although  we  have  only 

used  a  portion  of  the  information  available.  The  curl  operators 

may  be  taken  under  the  integral  sign,  because  the  field  point 
need  never  approach  the  fiber  surface. 


We  suppose  the  magnetic  field  on  the  surface  to  be  purely 


azimuthal,  so  that 


H+  (  r  '  ) 


0'  H ( s  '  ) 


(48) 


n'  x  H+  ( r ' )  =  s'  H( s  ’  )  . 


Using  the  identity  curl  curl  =  grad  div  -  div  grad,  the  first 
integrand  of  Eq.  (47)  takes  the  form 

s  •  Vx  Vx  n'  *  §.+  (r.')  g(kR)  =  s  *  (TV-  -V*V)s'  H(s')  g  ( kR ) 

=  H(s')  s  • VV-  s'  g ( kR )  +  k2(s-s')  H( s ' )  g ( kR ) 

=  H(s')  (s-D  (s'-v)  g  ( kR )  +  k2  (s  •  s')  H  (s')  g(kR) 

=  -H(s')  [a2/asas-  -  k2  (s  •  s')  g ( kR )  ,  (49) 


where  in  the  last  step  we  have  used  the  formal  notation  s  •  ^ 
d  I  d  s  and  s' '  •  V  =  ~  d  /  d  s  '  for  directional  derivatives.  Th  e  minus 
sign  arises  in  the  latter  case  because  the  primary  variable  has 
the  form 


R  =  _r(  s  )  -  r  (  s') 

For  the  second  integral  of  Eq.  (47),  the  electric  field  on 
the  surface  is  assumed  to  be  purely  axial,  i.e. 


s  •  VI  n'  X  E  +  ( r ' )  g(kR)  =  -s  •  (n'  x  E+)  x  Vg 

=  ( n '  x  E+)  •  ( s  x  vg ) 

=  E(s’)  (n'  x  s'*)  •  (s'  i  yg) 

=  -E(s')  (s  •  s'1)  (n*  •  Vg) 

=  -E(s*  )  <*  .  s'  )  ( n '  •  R)  k.  g '  , 

=  -E(s')  (s  •  s')  ( a / R) kg '  .  (51) 

Here  in  the  fourth  step  we  used  the  identity 

( a_  x  b_)  •  ( c_  x  d )  =  ( a_  •  c_)  ( b.  •  d )  -  (a.*  d )  (b_*  £) 

then  noted  that  the  scalar  product  n'  •  s'  vanishes  identically 
under  the  azimuthal  portion  of  the  surface  integration.  In  the 
last  step  above  g * ( kR )  =  dg/d(kR)  and  we  have  assumed  that 

A  A 

n'  •  R  =  a/R  .  (52) 


Note  from  Fig.  20  that  this  equality  only  holds  when  ls-s'l  <<  p, 
i.e.,  for  points  sufficiently  close  together  along  the  axis  that 


the  curvature  of  the  fiber  has  not  come  into  play.  No  approxi¬ 
mation  is  involved  here,  however,  because  the  integral  is  entirely 
negligible  otherwise,  as  we  will  now  see. 

Using  Eq.  (51),  the  second  term  of  Eq.  (47)  becomes 


.2  it 


+  b 


- ( k2 / 4 n  ) 


a  d0 


ds'  E(s'Hs-s')  (a/R)  g'(kR) 


-b 


.+  b 


=  - ( k/2 )  (ka)2  J  ds '  E(s ’ ) (s 


A 


s')  (1/kR)  g'(kR) 


-b 


+  CD 

r 


-  (l/2)a2  E ( s ) 


ds'[(s-s’)2  +  a2]~3/2 


E  (  s ) 


(  53  ) 


From  the  first  step  one  notes  that  because  ka  <<  1  only  that 

portion  of  the  integral  of  order  (l/ka)2  will  contribute  to  the 
final  result.  Noting  that 

(1/kR)  g'(kR)  -  - ( 1/ kR ) 3  =  -(1/k )3  [(s-s')2  +  &2]-3/2  , 


for  kR  <<  1,  the  remaining  steps  of  Eq .  (53)  are  straightforward. 
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Putting  the  results  of  Eqs.  (49)  and  (53)  back  in  Eq.  (47) 


and  carrying  out  the  azimuthal  integration,  one  finds  that 


+  b 


(  ia/2) 


ds ' H( s ' ) [02/as  3s'-k2(s. s' >Jg (s, s' )  -  E(s> 


-b 


=  -s  -Eioc  (s)  .  (54) 

We  can  express  this  result  as  an  integral  equation  for  total  line 
current  I(s)  by  writing 

I(s)  =  2  n  a i  H(s)  .  (55) 

and  introducing  a  surface  impedance  per  unit  length  given  by4 


Z  =  -(«n/2ff*kf)  J0  (kf a) /J0' (kf a)  .  (56) 

Here  kf  =  (w^jie  +  iupc)  1/2  is  the  complex  propagation  constant 
within  the  fiber,  JQ  is  the  Bessel  function  of  the  first  kind, 
and  we  have  assumed  that  both  E  and  H  within  the  fiber  vary  much 
more  rapidly  in  the  radial  than  the  axial  direction.  This 
assumption  is  consistent  with  the  requirement  of  moderate  to 
large  conductivity.  Note,  however,  that  if  one  were  to  represent 


I ( s)  as  a  Fourier  expansion,  then  eventually  the  axial  variations 


of  such  terms  would  dominate,  with  the  result  that  kf  in  Eq.  (56) 
would  have  to  be  modified. 4  That  is,  the  surface  impedance 
becomes  dependent  on  the  rate  of  axial  variation  of  current  when 
that  rate  is  large. 


The  tangential  E  field  can  now  be  expressed  as 


E  (  s  )  =  Z I  (  s  )  » 


I 


(  57  ) 


and  using  this  result,  along  with  Eq.  (55),  one  finally  obtains 


+  b 

r 


( 1/  4  Jt ) 


ds'  i(s')  ta2/as  as* 


-b 


k2( s. s' ) ]  g( s, s ' )  -  ZI( s) 


=  -s  •  Einc  (s) 


(58) 


Note  that  for  straight  fibers  s*s'  =1  and  this  equation  reduces 
to  the  usual  formula. 3  Also,  for  curved,  perfectly  conducting 
wires  Z - "0  and  one  again  finds  the  accepted  formula. 15 


4.3  Approximation  for  Special  Fibers 

For  the  general  case,  as  described  by  the  integral  Eq.  (58), 


it  is  clear  that  detailed  numerical  computations  are  required  in 
order  to  obtain  any  explicit  results.  If  the  fiber  axis  has 


radius  of  curvature  large  compared  to  wavelength,  or  is  made  up 


of  a  zigzag  series  of  straight  line  segments,  however,  then  the 
absorption  c r o s s - s e c t i on  is  readily  approximated  using  earlier 
results. 

First,  consider  the  case  where  the  fiber  axis  curves  only 
slowly  and  is  relatively  long,  compared  with  wavelength.  Then  we 
can  use  the  infinite  cylinder  result  to  obtain  the  absorption 
cross-section  per  unit  length  ca/£,  where  this  ratio  is  given 
by  Eq.  (30).  Note  that  ca  will  be  a  function  of  position  s  along 
the  curved  fiber,  in  that  oa  depends  on  the  angle  0O  formed  by 
the  incident  E  vector  and  the  local  tangent  to  the  fiber  axis. 
The  total  absorption  cross-section  Sabs  is  then  given  by 

ds  Ua  (*)/£]  .  (59) 


It  is  now  straightforward  to  obtain  results  for  toroidal  or 
C-shaped  fibers,  or  other  configurations  meeting  the  slowly- 
curving  limitation,  by  numerical  integration  of  Eq.  (59). 

For  a  fiber  made  up  of  zigzag  line  segments,  it  is  more 
appropriate  to  emply  the  quasi-static  approximation  to  the 
absorption  cross-section  oajjS  as  given  by  Eq .  (18).  We  have 
already  seen  by  comparison  with  other  computations  that  this 
formula  is  quite  accurate  for  all  fiber  lengths  provided  the 
modified  depolarization  factor  L'  of  Eq .  (29)  is  used  for  k0f >  1. 


Now  for  a  fiber  of  N  segments,  having  individual  lengths  , 
n  =  1,  2  ,  ...  ,  N,  the  total  absorption  cross-section  bs  i * 
given  s  imply  by 

r-<  N 

^-abs  =  2^.  °abs  sin2  Gn  '  60) 

n=  1 

where  the  factor  sin^0n  (angle  between  incident  E  vector  and  nth 
segment  axis)  is  included  because  the  original  Eq.  (18)  was 
specifically  for  broadside  incidence. 


It  is  useful  to 

note  tha  t  for 

both  Eq  s . 

(59)  and  (  6 ) 

the 

c  r  o  s  s  -  s  e  c  t  i  o  n  Sabs 

is  a  linear 

sum  of  the 

cross-sections 

of 

individual  segments.  Because  each  such  segment  behaves  precisely 
as  a  straight  fiber  ander  orientation  averaging,  one  concludes 
that  the  or ie nt a t i on- av e r a ge d  c r o s s- se c t i on 5Zab s  for  the  curved 
fiber  will  be  just  equal  to  that  of  a  straight  fiber  of  the  same 
total  length. 

5.0  TARGET  OBSCURATION 

In  this  section,  taking  into  consideration  both  the  diffuse 
cloud  scattering  properties  and  the  beam  extinction  properties, 
one  requires  the  optimum  conditions  for  target  obscuration  from  a 
mass  efficiency  standpoint,  over  a  wavelength  range  where  the 
theory  is  believed  to  be  most  accurate.  It  is  appropriate  to 
employ  the  quasistatic  approximation,  in  which  particles  are 
assumed  to  be  small  in  comparison  with  sensor  wavelength. 
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Consider  the  situation  when  a  target  of  (radar  or  optical) 
cross-section  aj  is  partially  concealed  by  a  cloud  of  particles, 
the  cloud  having  an  incoherent  cross-sectionSc.  Because  of  this 
incoherence,  the  cross-sections  are  additive  and  one  has 

Observed  cross-section  =  Sc  +  aj  e~^Vr  ,  (61) 

where  y  and  x  are  the  decay  constant  and  thickness  of  the  cloud, 
respectively. 

We  now  assume  that  the  sensor  cannot  detect  the  target  if 

£c  1  KoT  e'2rx  (62) 


where  K  is  a  figure  of  merit,  e.g,  for  K  =  0.1  the  signal  from 
the  target  would  be  10  dB  down  in  the  "noise"  of  the  cloud 
return.  The  diffuse  scattering  from  the  cloud  has  earlier  been 
computed  to  be5,6 

where  Vp  is  single-particle  volume,  the  reduced  conductivity  is 
given  by  ac'  =  k0e"  =  y/n0  /  e  0  oc  (in  terms  of  conductivity,  or  the 
imaginary  part  of  the  relative  dielectric  constant),  and  Ac  is 
the  geometrical  cross-section  of  the  cloud,  or  the  sensor  beam 
cross-section  at  the  cloud,  if  the  latter  should  be  smaller. 
Using  this  result,  and  taking  the  equality  in  Eq.  (62)  then  gives 


jf  \—J  VP  °c  =  K  0T  e-^TT  .  (64) 

The  decay  constant  for  the  cloud  is  given  by 


y  =  n  5ext  ,  (65) 


in  terms  of  the  number  density  of  particles  and  the  orientation- 
averaged  extinction  cross-section  per  particle.  Assuming  that 
absorption  effects  dominate,  one  has 


°  e xt  ~  °abs  * 


o  V 
c  p 


3 [ 1  +  (Le")23 


*  0C'  Vp/3 


(66) 


where  in  the  next-to-last  step  the  quasistatic  approximation  was 
employed,  and  in  the  last  step  we  noted  that  optimum  absorption 
will  occur  when 


Le"  =  4(a/£)2  [  n(£/a)  -1]  e"  <<  1  (67) 

(L  is  the  depolarizing  factor).  Note  also  that  the  total  mass  M 
of  particles  can  be  written 

M  =  p  Vp  5t  Ac  .  (68) 

Employing  Eq s .  (65),  (66)  and  (68)  in  Eq.  (64),  and  taking  the 


logarithm  of  the  result  determines  the  required  mass  to  be 


M 


3  p  A 
_ c 

2oc  ' 


in 


8  VT~  n 

/!li) 

V  l  V  n  9 

Lo  v  p  °  c 

\  Ac/ 

(69) 


The  quantities  K,  o-p  and  Ac  are  prescribed  by  the  logistics  of 
the  task  at  hand.  Thus,  to  minimize  M  one  first  must  seek  a 
material  with  smallest  possible  value  of  the  ratio  p/erc’.  In 
addition,  it  is  desirable  to  obtain  as  large  a  value  of  the 
product  k02  Vp  ac '  as  possible,  in  order  to  minimize  the  loga¬ 
rithmic  term  in  Eq.  (69). 


For  concreteness,  suppose  that 


Le"  =  0.1  . 


k0£  -  0.1  .  (70) 


From  Eq.  (67)  we  then  find  that 

0.01  ( £ / a  ) 

a  4  ac  1  Un(iU)  -1] 


(71) 


This  equation  can  be  used  to  determine  particle  radius,  once  the 
length  and  conductivity  are  known.  Typical  particle  design 


parameters  can  now  be  listed: 


r 


3 .x  10-2 


m 


C  =  4.8  x  10-4  m 
cc  =  5  xlO^  mho/m 
a  =  2 .2  x  10-8  n 

b )  infrared  region 
Y  *  3  x  10_5  m 
C  =  4.8  x  10-7  m 
oc  =  1  x  10^  mho/m 
a  =  1.0  x  10-8 

6.0  TAILORING  OF  PARTICLE  PARAMETERS  FOR  SPECIFIC  APPLICATIONS 

It  is  the  purpose  of  this  section  to  demonstrate  the  tech¬ 
niques  by  which  the  parameters  (length,  radius,  and  conductivity) 
of  thin  conductive  fibers  can  be  adjusted  so  that  a  cloud  of 
these  particles  will  have  selected  specified  spectral  charac¬ 
teristics.  These  characteristics  are:  (1)  large  absorption  and 

small  scattering,  (2)  large  scattering  and  low  absorption,  and 
(3)  transparency.  Two  frequencies  were  arbitrarily  chosen  for 
the  cases  to  be  analyzed.  These  are  f^  =  lO^-®  Hz  and  f2  =  10^ 
Hz,  corresponding  to  wavelengths  of  3  cm  (microwave)  and  300 
microns  (infrared). 

The  analyses  to  be  discussed  will  be  based  on  our  "extended 
quasistatic  model",  the  basis  of  which  is  described  in  Section 


2.2. 


As  demonstrated  in  Section  3,  this  analytical  model  produces 


results  which  are  in  surprisingly  good  agreement  with  the  results 
obtained  using  the  much  more  rigorous  variational  technique  of 
Section  2.1.  To  be  more  explicit,  the  quasistatic  model  can  be 
used  to  calculate  the  absorption  cross  section  over  a  very  wide 
range  of  k0C  ,  including  k0C  >>  1.  The  scattering  cross  section 
computations  resulting  from  the  quasistatic  model  are  valid  only 
in  the  range  k0C  jt . 

In  all  of  the  cases  analyzed,  the  results  of  the  quasistatic 
calculations  are  directly  compared  with  the  corresponding  results 
of  the  variational  technique.  We  find  these  comparisons  to  be 
quite  remarkable. 

6.1  Constitutive  Equations 

In  addition  to  the  equations  of  Sections  2.2,  2.4  and  2.5, 
the  following  equations  are  utilized  in  the  foregoing  analyses: 

Or i e nt a t i o n- av e r a g e d  absorption  cross-section 
(see  Eq.  21): 

ffabs  =  (3)  aabs  (72) 

0 r  i  e n t a t  i  o n- a v e r a g e d  scattering  cross-section  (see 
Eq.  22) : 

^sca  =  (5)  °sca  (73) 

High  frequency  depolarizing  factor  (see  Eq .  29): 


L'  =  4  I  fn  -1  ]  (k0£  1  n )  (74) 

The  inclusion  of  the  above  equation  permits  extension  of  the 
calculation  of  absorption  cross  section  well  into  the  k0t  >>  1 

regime.  Due  to  this,  we  refer  to  the  present  theory  as  the 
"extended  quasistatic  theory." 

In  this  treatment,  we  are  using  the  Drude  (free  electron) 
model  for  the  complex  dielectric  constant  (Section  2.4)  and  the 
Fuchs  model  for  the  dependence  of  conductivity  on  particle 
dimensions,  and  we  have  purposely  chosen  frequencies  and  u>2  such 
that  possible  anomalies  due  to  these  do  not  appear.  Also,  for 
purposes  of  analysis,  we  can  make  the  simplification  e  =  ie”. 
The  computer  program,  however,  does  not  utilize  this  simpli- 
f  ication. 

6.2  Absorb  at  3  cm.  Reflect  at  300  microns 

These  criteria  lead  to  the  following  mathematical  statements: 


At  frequency  idj  =  2n  x  10^-®  Hz 


(1) 

(2) 


<<  1 


(75  ) 


Let 


g  sp  ( 1 ) 
g  sp  ( 2  ) 


.01 


(76) 


(Le"  (1) )2  <<  l 


(77) 


Equation  (7)  ensures  that  the  absorption  cross  section  will 


be  at  its  maximum  value  at  frequency  cuj  ,  but  will  be  reduced  by  a 
large  factor  at  frequency  u>2  =  2  n  x  1012  Ilz  .  Also,  we  know  that 


the 

scattering  cross 

section  will 

be  at  least  as 

great 

a  t  u>2  a  s 

at  a>  i  .  This  will  cause  o 

s  ca ( ^ ) /oabs (2)  to  be 

large. 

which  is 

what 

we  seek  at  io2  . 

From  Eqs.  (76), 

(78)  , 

(72)  , 

and  (73),  we  obtain 

(with  the 

appropriate  definitiv 

e  equations 

of  Section  2.2) 

equations  for 

the 

conductivity  ac. 

radius 

a  and 

length  £  : 

e  c  2  2/3 

eoC  (.03) 

(79) 

u  c 

a2 

(10  [£n 

If)  - 

173  ' 

l]} 

a  ~  C  ( 

.03 ) 1/3 

,  .1/2 

<6(1 _  .  _  _ 

(80) 

,  \  1/2 

(10  ( £n 

(f)  - 

,..1/6  1/2 

11)  (oc) 

£  =  2a  jlO 

(f) 

-1] 

1/2 

/  "c  \) 

(  81  ) 

\  wi  l0lj 

when  Eq s .  (79)  and  (80)  are  equivalent,  i.e.,  we  can  choose  a 
conductivity  and  solve  for  a,  or  choose  a  radius,  a,  and  solve 
for  the  conductivity  ac.  We  will  do  the  former.  Knowing  both  a 
aod  °c»  we  then  calculate  the  length  from  (81).  Note  that  the 
term  in  square  brackets  in  (79),  (80),  and  (81)  is  very  insensitive 
to  large  variations  in  ( £ / a ) .  This  is  especially  true  when  it  is 
raised  to  fractional  powers  in  (79)  and  (80).  Therefore,  we  can 
come  quite  close  to  the  desired  results  if  we  simply  let  [  ]  =  2.5. 
This  permits  direct  estimation  of  a  and  £,  given  a  preselected 


value  of  a  r  , 


Table  I 


Co  n  d 

= 

1 . 00E+04 

Radius 

= 

8 . 6  8E-06 

Length 

- 

3 .7  IE-03 

k0  (  1 )  £/2 

= 

3  .  88E-01 

Cond 

1 . 0  0E+0  5 

Radius 

= 

2  . 7  5E-06 

Length 

3 .7  IE-03 

k0(l) C/2 

= 

3  .  88E-01 

Cond 

= 

1 . 00E+06 

Radius 

= 

8 . 6  8E-07 

Length 

= 

3 . 71E-03 

k0(l)  C/2 

= 

3 . 88E-01 

Cond 

1 . 00E+07 

Radius 

= 

2.75  E-07 

Length 

= 

3 . 7  IE-03 

k0  (  1 )  C/2 

= 

3  .  88E-01 

Cond 

= 

1 .00E+08 

Radius 

= 

8 . 6  8E-08 

Length 

= 

3  . 7  IE-03 

k0 ( 1)  C/2 

= 

3  .  88E-01 

Cond 

= 

1 .00E+09 

Radius 

= 

2  .  7  5  E-0  8 

Length 

= 

3  . 7  IE-03 

k0  ( 1 )  C/2 

= 

3  .  88E-01 

The  above  procedure  leads  to  sets  of  permissible  parameter 
values,  as  exemplified  in  Table  I.  Substitution  of  the  various 
parameter  sets  into  the  computer  program  yields  graphs  such  as 
shown  in  Fig.  21.  Since  these  graphs  are  identical  over  the 
wavelength  range  300  microns  to  3  cm,  only  one  is  shown.  This 
figure  is  representative  of  a  fairly  highly  conducting  metal 
(bulk  conductivity  =  10^  mho/m).  The  relaxation  time  x  is  that  of 
copper. 

From  Fig.  21,  we  see  that  (1)  indeed,  the  absorption  cross 
section  peaks  at  very  nearly  the  wavelength  =  3  cm,  (2)  that, 
in  the  case  of  absorption,  the  extended  quasistatic  theory  is  in 
excellent  agreement  with  the  variational  theory  for  wavelengths 
greater  than  about  50  microns,  and  (3)  the  scattering  as  calculated 
from  both  theories  is  in  excellent  agreement  when  k0f £  1. 


We  farther  conclude  that  the  mathematical  "design"  procedure 
demonstrated  in  this  sub-section,  although  not  analytically 
exact,  provides  good  parametric  sets  for  the  solution  of  the 
stated  problem,  since  the  objectives  are  met  at  the  two  specified 
wavelengths. 


6.3  Absorb  at  3  cm.  Absorb  at  300  microns 

We  have  quite  different  criteria  for  this  problem  as  compared 
with  the  preceding  problem.  Here,  we  wish  to  have  a  very  highly 
absorbing  cloud  which  has  low  scattering  cross  section. 


In  this  case,  the  cloud  scattering  cross  section 


given  by 


(82) 


where  Ac  =  projected  area  of  the  cloud  or  that  portion  thereof 
which  occupies  the  solid  angle  of  the  incident  beam. 

Obviously,  we  wish  to  make  the  absorption  cross  section  of 
the  particle  much  larger  than  its  scattering  cross  section.  This 
must  be  true  over  at  least  the  wavelength  range  from  300  microns 
to  3  cm. 

The  above  considerations  lead  to  the  following  mathematical 
criteria  : 


(1)  In  order  to  have  maximum  absorption  over  the  specified 
wavelength  range,  we  want  Le"(l)  <<  1  and  Le"(2)  <<  1. 


This  means  that  the  particles  should  be  very  thin  and 
(perhaps)  not  too  highly  conducting. 

(2)  Since,  under  the  above  conditions,  5sca  ~  oc^  and 

cabs  ~  Vp  oc,  the  ratio  (osca/°ab s)  is  proportional  to 
k0^  Vp  oc.  Therefore,  in  order  to  maintain  high  cloud 
absorption  and  low  cloud  scattering,  we  wish  to  keep 
the  product  Vp  oc  low  (we  have  no  control  over  kD ) . 

A  little  bit  of  experimenting  with  the  extended  quasistatic 
( EQ )  computer  program  yields  appropriate  sets  of  parameters.  The 
parameter  set  corresponding  to  minimal  total  mass  is  a  =  100 

Angstroms,  £  =  100  microns,  and  erc  =  10^  mho  /  m  as  the  bulk  conduc¬ 
tivity.  The  results  of  using  these  parameters  in  the  EQ  and 
variational  codes  are  shown  in  Fig.  22.  One  can  trade  off  a 
larger  radius  for  a  lower  conductivity  and  achieve  similar 
results,  but  with  a  somewhat  higher  ratio  of  (B sca /oabs )  • 

Note  again  the  remarkable  agreement  between  the  computations 
based  on  the  two  independent  theories! 

6.4  Reflecting  Modes 

It  is  easy  to  design  particles  having  very  low  absorption 
and  relatively  high  scattering.  Experience  has  shown  us  that  if, 
for  example,  we  choose  a  highly  conducting  metal  such  as  copper 
or  aluminum,  and  select  a  radius  of  one  micron  or  greater,  the 
particle  will  be  essentially  completely  reflecting  for  at  least 
all  wavelengths  below  30  microns.  le  can  tailor  the  scattering 


cross  section  to  become  proportional  to  u1 2 3 4  at  wavelengths  signif¬ 
icantly  below  that  for  which  kD  ~  1,  which  of  coarse,  represents 
Rayleigh  scattering  from  "perfectly  conducting"  wires.  Thus, 
the  aggregate  of  particles  can  easily  be  made  to  be  reflecting  at 
one  wavelength  and  essentially  transparent  at  some  (significantly) 
lower  wavelength.  Chaff  clouds  behave  in  this  way. 

6.5  Transparent  at  3  cm.  Absorb  at  300  Microns 

In  this  case,  we  want  the  scattering  cross  section  to  be 
much  lower  than  the  absorption  cross  section  at  all  wavelengths 
under  consideration.  The  requirements  due  to  this  and  the 
transparency  criterion  at  <d^  are  given  below: 

(1)  Make  Vp  ac  as  small  as  possible,  consistent  with  the 
other  criteria. 

(2)  Make  Le  "  =  1  at  u>2  . 

(3)  Make  kQX  <  1  at  u>2  • 

Application  of  these  criteria  leads  to  the  parameters  listed 
in  Figure  23.  Note  that,  at  k  =  3  cm,  the  absorption  is  three 
orders  of  magnitude  lower  than  that  a t  X  =  300  microns.  The 
scattering  throughout  the  <o2  to  ui2  spectrum  is  much  lower  than 
the  absorption. 


Comments  Relative  to  Particle  Parameter  Tailoring 

(1)  Oar  extended  quasistatic  treatment  provides  a  very 


useful  tool  in  the  design  of  particles  for  specified 
applications.  In  some  instances,  the  e.q.s.  theory 
does  not  yield  accurate  quantitative  data  and  should  be 
used  primarily  as  a  first  step  to  be  followed  by  the 
full  computation  using  our  Variational  method. 

(2)  The  calculation  of  scattering  cross  sections  using  the 

EQ  theory  is  limited  to  the  range  k0f  /2  <.  1.  Unfor¬ 

tunately,  the  accuracy  of  the  scattering  cross  sections 
using  the  variational  theory  decreases  substantially 
for  k0i  >>  1,  This  is  discussed  in  Section  3.  Our 
numerical  integration  technique  is  presently  undergoing 
substantial  improvement.  Also,  we  are  presently  using 
the  Hewlett  Packard  HP  9000,  Model  520  computer  for 
these  calculations  -  providing  much  higher  computation 
spe  e  d. * 

(3)  We  have  included  backscatter  cross  section  computations 

in  Figures  21,  22  and  23  for  those  applications  that 

require  this. 

(4)  Within  the  bounds  of  present  technology,  it  is  possible 
to  create  particle  parameters  that  could  provide  useful 
spectral  characteristics  over  a  wide  wavelength  range. 


*Note  -  added  in  proof:  At  the  time  of  submission  of  this 

report,  we  have  made  substantial  improvements  in  our  numerical 
integration  routines.  This  includes  k0  dependent  incrementing 
of  the  number  of  the  scattering  cross  section  computations.  The 
minimum  incrementation  of  scattering  angle  is  0.05  deg.  and  the 
maximum  is  1.0  deg.  This  gives  good  results  for  k0i  i.  350. 
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